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I. INTRODUCTION 



A. Closure theories and formalisms in turbulence 



In turbulence, statistical closures are techniques employed to close the moment-hierarchy established in a statistical 
treatment of the Navier-Stokes equation (NSE). Such theories postulate a relationship between high- and low-order 
moments by way of physical arguments. The goal is to accurately describe and predict the statistics of a turbulent 
system while maintaining a strong connection to the underlying dynamics of the NSE. Closure theories can, in 
principle, allow efficient computation of turbulent statistics without the computationally intensive demands of Direct 
Numerical Simulation (DNS) This feature makes closure theories particularly attractive when computational 
power is limited. Although the ability to compute the full Navier-Stokes equations using DNS is currently and 
increasingly more tractable, closure-based computations are still able to provide useful insights into turbulent systems 
at a much smaller computational cost. 

The analytic study of turbulence can be classified into different formalisms from which specific closure theories can be 
constructed. Two key formalisms which have developed, and are the focus of this paper, are the Wyld and Martin, 
Siggia, and Rose (MSR) 0] formalisms. The closure theories themselves that can be obtained from these formalisms 
can be further classified into two main groups, Eulerian and Lagrangian, distinguished by the reference frame from 
which the dynamics is described. The former describes a fluid from a 'lab', where the fluid moves relative to a fixed 
frame of reference outside the fluid. In particular this is the frame of the Navier-Stokes equations. The Lagrangian 
description is in fact a re-formulation of fluid dynamics that calculates statistics by following fluid particles. The 
work presented here docs not consider the Lagrangian closures but instead focuses on those of the standard Eulerian 
formulation. 

One of the earliest closures for homogeneous isotropic turbulence was the quasi-normal approximation. In this 
approximation the fourth-order correlation is written in terms of products of second-order correlations 0, HI . However 
turbulence has highly non-Gaussian correlations and so the resultant predictions from this simple approximation were 
inadequate, in fact leading to total kinetic energy having negative values The failure of quasi-normality triggered 
research that made improvements upon it, resulting in the EDQNM-mcthods (sec for example Q), which are a subset 
of the more general class of Renormalized Perturbation Theories (RPTs). The RPT approach re-sums certain selected 
infinite terms from the conventional perturbation expansion. The underlying idea is these resummations can capture 
essential non-perturbative physics of turbulence. There is a body of literature that has grown up around this work 
and some principle sources elaborating these directions are the books by Leslie [1] and McComb [1]. 



B. Equations and quantities 

The main focus of the theories described above is to compute the two-point, two-time correlation function, 
(uq(x, t)7i^(x', i')), of a homogeneous and isotropic field of incompressible fluid turbulence. Although this quantity 
is defined in real-space, we will find it easier to start with the spectral Navier-Stokes equation (NSE), 

j-l-l=k 

defined using the following Fourier convention 

a„(x,t) = /"dka„(k,t)e""^, a„(k,t) = -i^ /'dxa„(x,Oe-"'", (2) 



where, for the velocity field a = u and for the forcing a = f. The notation j -I- 1 = k under the integral signs indicates 
that the integration variables are constrained under convolution. Note that the pressure term has been removed using 
the continuity equation, k • u(k, t) ~ and that this results in the introduction of the tensor M^°^^(k). Formally, it 
is called the 'momentum transfer operator' (see and is defined as 

M^:M ^ ^ (k,S^, - ^ + k,<5„, - ^p-) (3) 



UkpP^,{k)+k,P^p{k)\ (4) 
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The tensors in the last hne are a result of the isotropy and are defined to be 

Pa,(k)=J„,-^. (5) 

We use this formulation of the NSE to then derive an expression for the spectral correlation function, 

C^^,{k;t,t')S{k + k') = {u^{k,t)u^,{k',t')), (6) 

which we can later inverse-transform to real-space. However, it is well-known that the nonlinear term prevents a 
statistical treatment of the NSE from giving a closed set of equations; an equation for the (two-point) correlation 
function requires knowledge of the three-point correlation function, and so on. The formalisms described in this paper 
are attempts to circumvent this problem in a systematic way. We first turn to a famous example of a Renormalized 
Perturbation Theory. 



C. Example: The Direct Interaction Approximation 



Kraichnan's Direct Interaction Approximation (DIA) |10l - |l2l | pioneered the renormalized perturbation theory ap- 
proach to homogeneous isotropic turbulence. The basic idea of this approach was to focus on two quantities, the 
velocity correlation function and the response propagator function, and obtain approximate equations to calculate 
them. 

The main hypothesis of the DIA was that the wave- vector interactions in the nonlinear term would be dominated 
by a single triad satisfying k = j + 1. Kraichnan exploited this concept and used it to create a unique perturbation 
expansion that could be used to bring about a closure to the statistical hierarchy. Ultimately, he derived a closed 
set of equations that use only these so-called directly interacting wave- vectors. In the final form, the set includes an 
equation of motion for the exact correlation function of velocity-field coefficients. 



(dt + lyk^) Caa' (k; t, t') = M^'>^{k) J d^i dH 



j+l=k 



2 / ds Rp0,{j; t, s)Ml,Vs,{j)C^s{\; t, s)C,„-(-k; s,t') 



ds Ra'd-k; t', s)M'^l,^, (k)C^y (1; i, s)Cp0, (j; s) | , (7) 

along with an evolution equation for a quantity known as the renormalized propagator function, i?Q,Q' (k;t, i'), 

(at + i/k2)i?,„,(k;t,i')-4M™^(k) jd^]dh X 

j+l=k 

Ids Rpp, (j; t, s)M<';;,3(j)C^y (1; t, (-k; .s, t') \ 



t' 



Paa'{k)5{t~t'). (8) 



We give these equations without derivation so that comparisons can be made later in the following sections: more 
information can be found in his original papers as well as in BeranfH, McCombdH, Kida and Gotoh jl5|, and 



Krommes [16|. Leslie's book |8| is largely dedicated to Kraichnan's works from this period and provides many insights. 

The DIA, although successful in low- to moderate- Reynolds numbers, fails to produce a Kolmogorov inertial range. 
Kraichnan himself showed that the DIA gives an inertial range with k~^/^ [ll| and suggested that the DIA does not 
properly deal with the indirect interactions [13 ' manifest by the DIA's inability to decouple the large-scales from 
the viscous scales BBH- As will later be seen the momentum transfer terms are in effect vertex functions. The 
indirect interactions are intrinsically associated with these vertex functions leading to the notion that "the whole 
problem of strong turbulence is contained in a proper treatment of vertex renormalization" Q . 

The success and failings of the DIA led to further closures based on renormalized perturbation theories. Notable 
ones are those of Nakano [l^ , and the Local Energy Transfer (LET) theory of McComb [1, [20l - [30| ; the latter being 
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the only purely Eulerian theory which is compatible with the Kolmogorov inertial range. Convinced of the perceived 
intrinsic failings of the DIA based on an Eulerian framework, Kraichnan reformulated fluid dynamics to use Lagr angian 
variables and produced the Lagrangian-DIA [sH ]. This also led to many off-shoots notably those of Kraichnan [sl. lol. l32l| 
and the LRA of Kaneda [s^ [S^] and Kida and Goto [l^ . The Eulerian-DIA still enjoys some use notably in the 
regularized-DIA (RDIA) of Frederiksen [ssj - issf . For completeness, we should mention the separate line of development, 
the functional formalism of Hopf [s^ , leading to the self-consistent theories of Edwards [iol and Herring [U \^ . 

The sections to follow detail the two formalisms of Wyld and MSR that aim to achieve a statistical closure which 
properly deals with both the direct and indirect interactions. Both formalisms can be applied more generally to 
classical dynamics systems, but in their original papers were applied to turbulence. It will be argued that both 
formalisms, although having different methodology, are equivalent; it will be explicitly shown that both formalisms 
obtain at lowest order the DIA and also agree up to, and including, fourth order in the pcrturbativc expansion. 



II. THE WYLD FORMALISM 

The W yld formalism Q is a perturbative analysis of the statistical turbulence. It represents one of the earlier 
attempts [l9l l43l| to extend the methods of quantum field theory (QFT), specifically those of diagrammatic represen- 
tation, to the problem of classical turbulence. 

In the approach used by Wyld, the velocity-field is expanded in a perturbation series, with an associated dia- 
grammatic representation, which then is used in a statistical average to obtain the two-point velocity correlation 
function. As the perturbation series are in fact infinite, a systematic renormalization method is employed to reduce 
the series into a more manageable form. This results in integral equations, which at lowest non-trivial order, repro- 
duce the Kraichnan DIA result. We show in subsect. Ill CI that Wyld's original renormalization procedure suffered 
double-counting issues as noted by Lee [4^ , and offer a formal procedure that can remedy this problem. 



A. Wyld's Perturbation Method 



The main focus in examining Wyld's method is the renormalization procedure. However, the fundamental or 'bare' 
equations must be established prior to renormalization. The following briefly explains Wyld's construction of the 
velocity correlation function via a perturbation expansion of the velocity field. There are a few places in which the 
following summary of Wyld's method deviates from his original work and this will be pointed out below. However, 
these differences do not change in any essential way his original analysis. 



1. Wyld's perturbation expansion 



Once again, we start with the NSE in Fourier space, 

(^dt + uk^)u^(k,t) - Uk,t)+AQ^{k) jcP] up{lt)u^{\^^lt). (9) 

It must be pointed out that our approach already differs somewhat from Wyld's in that he also Fourier transformed 
the time variable. As such, in the Wyld analysis the wave-vector k and wave-frequency to are then lumped together 
into a 4- vector k = (k, — w) and the tensorial NSE is abandoned in favor of a simpler one-dimensional 'model' equation. 
We will not make any of these changes, and rather work directly with the 3d NSE. However we will follow the basic 
formalism set up by Wyld. 

Inverting the linear differential operator on the LHS of © to the RHS results iir the following form of the NSE, 

t 

— oo 

t 

— OO 

(10) 



5 



A bookkeeping parameter, A, has been multiplied to the nonlinear term for the purposes of the perturbation expansion; 
it will later be set equal to unity. The next step is to consider a perturbation expansion of the NSE, 

(k, t) = u^°\k, t) + Xn'^\k, t) + A2u^^)(k, t) + . . . , (11) 

This can be substituted for each velocity field in Eq. (|10p and then expressions can be matched by powers of A. At 
the lowest order, 



t 



<>(k,i) = y'dt'<i'(k; t,i')/c'(k,t'), (12) 

and one can establish a response function i?^"^(k; t') such that 

(dt + lyk^) [i?™(k; t,t')] - Pa/3(k)J(i - t'). (13) 

Already, it may be seen that there are many variables, arguments, and indices to keep track of, therefore it is useful 
here to introduce a reduced notation: 

u^^\k,t)-^ul^ (14a) 
i?™(k;f,t')^i?r (14b) 

<^^(k)/d'j -^M^. (14c) 

This notation will be less cumbersome for the reader to follow; vector indices and time arguments can be determined 
later where needed. The spectral NSE in the new notation becomes 

{dt + iyk')uk = h + XM^"'uiUk-s, (15) 
and the equations, Eq. ((TO)) and pT|) . used for the perturbative treatment are now respectively. 



uk = i?L7k + A<'M^ '"j^k-j (16) 



~(0) . \ ~(1) , -i2~(2) . \3~(3) /1'7^ 

Uk = "k + ^"k + X U^^ + \ U^^' . . . (17) 



In using this notation, the integral following the momentum transfer operator is always a convolution, where the wave- 
vector arguments of the convoluted functions must add up to the wave-vector of the momentum transfer operator 
immediately preceding them. Some care may be initially needed to keep track of these integrated wave- vectors; a 
simple rule that adjusts for this is that all non-k wave-vectors are integrated out. 
The perturbation terms by order in A are 

A° : = (18a) 

A^: w'^' =A<'Mj'"w^'w<°'_j (18b) 

A^: tf^>=A2<'Ar(2^fS<^'=) (18c) 



The term w'^Lj in the integrand of the expression for w'^' can be replaced by its definition, leaving u*^' written only 
in terms of u'-°\ In fact, any order i?"' may be written as product of (n + 1) fl^^'s. For example, the last term in the 
above equation for the perturbation expansion of u*^' can now be written as 

A^ : = 2X^R'^'M^">R'°l.M^"2iufu'°'ull._,. (19) 
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2. Wyld's correlation 

Wyld approaches the correlation of two velocity fields very simply by considering the average of the product of two 
velocity field expansions of u's: 

(?i„(k, t)7i^(k', t')) = ^ {<\K t) + A^iL"(k, t) + • • • ) X 

[u^°\k',t') + Xul^\k',t')+ ■■■)y (20) 

The zeroth-order and exact correlators are given respectively by 

cz{\^;t,mk'+k) EE (unk,t)ui;"(k',i')) = K'KD, (21) 

C^Jk;t,t')S{k' + k) EE {u^{k,t)ujk',t')) = {u^u^,), (22) 
and their reduced-notation counterparts are 

C™(k;i,i')W ^C^rW, (23) 

C„„(k;i,OW -^C^k4+k'- (24) 

The delta-function on the LHS of these definitions is a result of the construction of the Fourier-transform of the 
real-space correlation equation, 

(7i„(k,t)^^(k',t')) = y^'x Jd^u^{^,t)uU^ + v,t')e-''''-^e-'''<^+'^^^ 

= (2^)6 J'^''' yd'r(^i„(O,0wc.(r,t'))e-^^'''+*'^ "e-"' - 
= (2^ jd^^CU^-t,t')e-^'^-^5{\^ + k) 

= C^^{k-Xt')6{k! + k). (25) 

The second line uses the homogeneity constraint, ({{^(x, t){2(^(x -|- r, i')) ~ ({{^(O, t)u^(r, t')). 

Note that the zeroth-order velocity field expansion terms are random Gaussian functions as is implied by the 
temporal delta- function in (|13p . Recall the correlation of an odd- numbered product of random Gaussian variables 
vanishes. Thus in our analysis this means all terms in the perturbation expansion with an odd power of A will vanish. 

Continuing in reduced notation (without odd-order moments), a series-expansion for the exact correlator is obtained, 

+ A^(("'k'«'k;) + ("'k'"'k;) + K'^k;)) 

+ A^(«^'k!) + ("'k'^'k^) + (^'k'^'kO + (^<k'«'k!) + ("'k'^'k^)) 

+ C'(A^). (26) 

As mentioned above, all terms can be written as products of zeroth-order velocity fields. For example consider the 
last of the second-order correlations (or moments): 

(«'k'0 = 2i?L°'M^°><ijM^!.'j(Gf <'u<°'_j_,«<^!) . (27) 

Another property of random-Gaussian variables is that any n-th order moment may be decomposed into a sum of 
products of lesser-order moments. In the above case, the fourth-order moment is decomposed into three pairs of 
second-order moments, 

(28) 
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Note that all possible combinations of second-order moments are created in this decomposition. What is immediately 
useful here is that the fourth-order moment can be written as pairs of second-order moments, which more importantly 
are zeroth-order correlation functions. 

Using the definition of the (zeroth-order) correlator, the above equation becomes 

The last term vanishes as it implies AI^°2j^k-}, which vanishes by definition. Cleaning up leaves 



A similar calculation can be made for the other terms, giving the correlation equation to second-order, 

1(0; 



Ckf^k+k' = C'k'''^k+k' 



+ 2R'^^ M^"' R'°^ M^V q"' C'^K 5k+k' 

+ 0{X^). (31) 

This can be applied to all orders. However, as this is an infinite expansion, a full calculation will be intractable. This 
led Wyld to use a diagrammatic resummation, an approach that would contain the effects of all orders generated by 
the perturbation expansions into a more manageable set of equations. 

B. Wyld's Diagrammatic Method 

Wyld associated diagrams to the terms of the perturbation expansion. These diagrams could then be combined 
to form order-by-order the velocity correlation. This procedure in turn produced a graphical expansion for the exact 
correlator function. In this subsection this diagrammatic formulation is presented and the resulting expressions for 
correlation function based on these diagrams is given. 

1. Defining diagrams 

Wyld diagrammatic notation assigns particular symbols to the various terms present in the perturbation expansion 
of the velocity field: 



t__k t' 
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O i?^"'(k; i,t') ^ i?<;" Bare propagator (33) 



^ ^Ca-ylk) ^ A^k°' Bare vertex. (34) 

T k-j ' 

These are placed into the relevant equations for the perturbed expressions. 

Following the example of Eq. (|19p . the second-order velocity term in the perturbation expansion, 

~(2) _ r,p(0) l\/riO) ~(0) p(0) ,^(0) f~m~m \ /or\ 

"k - ^k -'^k-j^k-j {^i "k-j-ij' (35) 



can be written in diagrams as 



-(0) 
u\ ' 

(0) 
k-j 

k-j-i 
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The effectiveness of tliis diagrammatic approach is appreciated by seeing it with all variables, arguments, and indices 
restored. 



i?™,(k;t,f 



(0) 



(k) 



J?'"',(k-j;t',t")' 

ft 



u<o>(k-j-l,t") 



for which the corresponding analytic equation is 



t 

— oo 

t' 

jdt"R^^l^,{\.-]- t',t")M^^U^-i) jdh (uT{\,t")uf\\.-i-U")). (36) 

— oo 



2. 'Tree' Diagrams 



Tree diagrams in the perturbation expansion always begin on a propagator and end in external velocity legs. The 
equations and their diagrammatic representations for these terms are given below to fourth-order: 



'k-j 



4°'. 



(37) 



u 



"k-j-1 



(38) 



X". 



g 



k-j-l-m 



(39) 
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"h "l-h"m"k-j-l-mj 



' ~(0) ~(0) ~(0) ~(0) ~(0) 



k-j k-j k-j-1 k-j-1 k-j-l-m k-j-l-m ' 



-(0) -(0) 

) 

i(0) 



\ -(0 



) 

k-j-I-r 



-(0) 

■ 8i?<°' Mi"' i?L°' M<°' i?<°' , M<°' i?<°i , „ M'°' X I 8 



i" ' ■"m^n' '"k-j-l-m-n ) 



) 

k-j-I-r 



-(0) 

"m-i™ • (40) 

Due to their apparent shape, these are the so-called "tree" diagrams, with single propagator lines on the left as trunks 
and propagators ending in zcroth order velocity lines as branches on the right. 



3. Correlation diagrams 

Correlation diagrams arise from attaching tree-level diagrams together by fusing the velocity field terms at the ends 
of branches. These become the zeroth-order correlation functions. The diagram for the zeroth-order correlation term 
or 'bare correlator' is given by 



O C™(k;i,t') -> C^. (41) 



A similar diagram is used for the exact correlation 

t k t' 



|3 



O C„^(k;i,<') ^ Ck, (42) 



where the double lines are used to distinguish it from its bare counterpart. 

To see how the diagrams operate, it is instructive to examine the construction of second-order correlation terms 
from the tree- level diagrams. The first term considered here is the last of the second-order terms in Eq. (pS)) . 

= 4<" M^°> B!^ q"' C^",' <Sk+k' . (43) 

Diagrammatically, this corresponds to 



(0) 



\ -(0 



,(0, 
-(0) 



(44) 



Note that we will treat many of the diagrams in this paper as equations. An extra factor of 2 arises from the 
combinatorics, analogous to the Wick contractions in QFT of the fourth-order moment into products of second-order 
moments; for example we consider the combinations of the second-order term from above 



\ 



"k-j-1 "k-j-1 "k-j-r (45) 
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which results in 

K'n'^,) ^ {ufu'nKl,^K'') + + ("r"k!)«"'k-j-i) • (46) 

The last term above vanishes as in Eq. ()29p , introducing a rule that any diagram with a closed loop that is connected 
to the diagram by a single propagator line will vanish. 

The other diagrams in the second-order terms are obtained by a similar construction, 

= 4R'°^ M^°' R'^Ly M^°ij, C^"' C<°'_j, (5k+k' (47) 

implies 

q(0) 

AO)/ 



and 



V-j'-i' (48) 

("k%') -^k Mc -'i^k'Mc' ("j "k-j"j' "k'-j'/ 

= 2i?^''' M^°' i?<^°; M^°' q"' C^' 5k+k' (49) 



implies 




This procedure can be applied to all orders, and it can be shown that a one-to-one correspondence is established 
between diagrams and their analytical counterparts, with the correct numerical prefactors. A reproduction of primitive 
correlator diagrams to fourth-order is given in Fig. [TJ 

In Wyld's perturbation method all correlations of bare quantities involve Gaussian statistics but at all orders in 
the perturbation expansions it is formally exact. However this means retaining an infinite number of terms in the 
expansion of the correlation function. The next subsection sees the systematic renormalization of these terms into a 
manageable formula. 
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FIG. 1. Wyld's diagrams representing the correlator expansion up to and including the fourth-order terms. The labels shown 
on each diagram correspond to those given in Figure 2 of the original paper Q] . 
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C. Improved Wyld-Lee Renormalized Perturbation Theory 



Wyld's method of renormalization is a resummation of diagrams based on the emergence and recurrence of funda- 
mental, irreducible diagram units. These are so-defined since they are diagrams of a given order which can be found 
as parts of diagrams at all higher orders. To illustrate this, we use diagram (W6) in Figure [T] as an example; this 
diagram is essentially diagram (W2) in which the top correlator has been replaced with diagram (W2); the result 
is a reducible diagram of 0{X^) composed of an irreducible diagram of 0{\^) embedded into itself. This will be 
applied not only to correlators but to the propagator and vertex functions as well, introducing additional renormal- 
ized equations for these quantities. By this process, the resulting system of equations is constructed only in terms 
of renormalized quantities, establishing a new renormalized perturbation expansion. An analogy can be made with 
finding the irreducible diagrams in other diagrammatic methods, most appropriately those in particle physics [45| . 

The main controversies associated with Wyld's formalism are due to the renormalization, and therefore, the pro- 
cedure here will be different from Wyld's. In Wyld's approach for the correlator, he directly computes the 2-point 
velocity correlator, inputting the perturbation expansion at each order for the two velocity fields. Instead, we will 
deduce the two-point velocity correlation function through the Navier-Stokes equation. In particular, the new starting 
point is that which was used by Kraichnan for DIA [l^ (and McComb for the Local Energy Transfer (LET) theory [^). 
Namely we look at the spectral NSE multiplied by a second velocity-field coefficient, Uq,'(— k, t') and then averaged, 



With the definition of the exact correlator, Eq. ([2^ . the above equation can be viewed as an evolution equation for 
the exact correlator, eventually leading to the spectral energy equation. To compute the above equation, we apply the 
perturbation technique to the terms on the RHS, the force- velocity and the three-point velocity correlation functions. 
This approach was not originally used by Kraichnan but was found as an alternate route to deriving the DIA by Leslie 



When this correlation expansion is renormalized, there are two ways it differs from Wyld. First when the linear 
operator on the LHS is taken to the RHS, it remains a bare response function (see on the LHS of Eq. ((FT|) in 
its inverted form) and this term within the renormalization procedure will always remain unrenormalized. Second, 
the vertex associated with the momentum-transfer operator M^p^Ck) is outside the average, and will also not be 
included in the resummation. By our above procedure both these points are accommodated in a natural way into 
the formulation. Part of this renormalization procedure has been previously pointed out in the book by McComb 0] , 
where the first of the above two steps was noted. Keeping the vertex function on the RHS of Eq. ((5T|) in its bare 
form, is being noted here for the first time. As we will discuss below, both these steps are necessary to avoid double 
counting effects in the renormalization procedure. 

The process of renormalization firstly requires the identification of reducible and irreducible diagrams. The first 
criterion by which to classify diagrams accordingly is the ability to separate a diagram into two parts by severing a 
single correlator. Diagrams that can be split into two separate diagrams by cutting a single correlator are labelled 
by Wyld as 'Class- A'. We introduce a modification to the procedure here which is to further distinguish two types of 
Class-A diagrams: those diagrams with the correlators on the LHS of the left-most vertex, and those with separable 
diagrams connected by a single bare correlator that occurs on the RHS of the left-most vertex. These are labelled 
Class- Ai and Class- Aj^, respectively. 

For purposes that will soon become apparent, we note that bare correlators are given analytically by 



which reveals the propagators within each. This identification allows the derivation of the propagator diagrams to 
come. An example of cutting an external correlator is given diagrammatically by 




(51) 



i 



««'(k,t)^L"'(k',i')) = «^'(k; t,s)U,{k,s)fl,{k\s')RZ,{k'; i',.')), 



(52) 




with its analytic counterpart given by 



("k ' ) Rr < 



) 
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This diagram, (W4), is a member of the Class- diagrams; note that the equation still reads left to right. Using 
Fig. 2.2 as a reference, the Class-A^ diagrams are: 4, 17-25(odd), 26, 29-33(odd); and the Class-A^ diagrams are: 
3, 5, 16-24(even),27, 28-32(even). While it is not evident, it must be pointed out that the set of Class-A^ diagrams 
is smaller than its complement but without this classification diagrams will be created redundantly as was found 
already at fourth-order by WyldQ and Lee It is for this reason that we have introduced this subdivision of 

Wyld's Class-A diagrams. 

Note that Class- A^ diagrams can be written with the leftmost correlator in the form of a force-force correlation. 
We introduce an equation for these diagrams to fourth-order, first graphically followed by its analytic counterpart. 



(/k/k')= 



^(/k./k')- 



4--(/k/k')-^- 



0(A4), 



(53) 



i?r(/k/k')i?k' = <'(/k/k')<' + ^K'ifMK'K'K'RsK'c^' + 



(54) 



In writing these equations, a new function, the exact propagator, has been introduced and a new diagram is associated 
with it. 



0(A4), 



(55) 



?(0) 

i^k 



Ry = R 

The complete expansion of the propagator to fourth-order is given below 



(56) 



-I- 4 



(Wp3) 

+ S — W^^-f-r.-.. + 16 — V^^'^.r-., + 16— vO^^ 



(Wp2) 



(Wp4) 



/ (Wp.">) 



/ (Wp6) 



/(Wp7) 



-I- 16- 



-I- 16--- 



+ 16- 



(Wpll) 



/ (Wp9) 



-I- 16 



(57) 



Note that while the method used here to obtain these terms is different from Wyld's, in the end the same terms 
are obtained for the expansion. In determining this expansion, we have only considered a subset of Wyld's Class-A 
diagrams and enforced that these diagrams have a bare propagator on their left, breaking the leftmost correlation 
function that characterize the Class-Ao diagrams. As suggested by McComb Q, this maintains consistency with the 
DIA derivation of Kraichnan [1, [Ullill. A similar prescription will be used in determining the renormalized equation 
for the propagator as we will later see. 

The procedure for dealing with Class-A/j diagrams will be postponed as it is similar to that used to re-sum the 
remaining diagrams. Diagrams not classified as Class-A are designated as Class-B diagrams; these are further classified 
into reducible and irreducible based on finding embedded elements of low-order within diagrams of higher-order. An 
example of this can be seen by examining the diagram W6 in Fig. [1] 



It is readily seen that the elements between the two outermost vertices is the diagram W2 in Fig. [T] this is also 
equivalent to the W2 diagram by replacing the top correlator with itself. This will be given in more detail below 
however it will be useful to include another function at this time. 
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Wyld introduced an exact vertex function as an expansion without giving a detailed account of how it was derived 
but instead wrote a diagram series expansion for the exact vertex function, 



4iX 



0(A5) 



(58) 



4M^°' C^"' R^' A'/J"' Mj^"' 



AM^°^ i?!"' R'°l M^°^ C'"' M^"' 



O(A^). 



The renormahzed vertex diagram is defined in the same manner as its bare counterpart, 



(59) 



^ k-j 



(60) 



Starting with these irreducible diagrams, Wyld generated the full expansion by replacing any of the constituent 
elements with a higher-order element. For example, replacing the correlator, a propagator, or a vertex in the first 
term of Eq. (|58|) results in the following terms, respectively: 



Irreducible diagrams can then be found by selectively removing correlators, propagators, and vertex corrections from 
Class-B diagrams. The set is left with diagrams that cannot be constructed from non-trivial (bare correlators, 
propagators, or vertices) elements. Under this classification, there are two up to fourth order. 




(W2) ■ 



(W7) 



Starting with the second-order Class-B irreducible diagrams, all Class-B diagrams (except the Class-B irreducible 
diagrams that arise at higher orders) can be generated by replacing the appropriate correlator, propagator, and vertex 
corrections with their respective expansions. This is demonstrated using three correlator diagrams. 



(W6) 



o- 




The three correlator diagrams are constructed as follows: 
a. Construction with a Correlator Correction 



-o- 



+ 



+ 



(61) 



Replacing a bare correlator with an exact one in the irreducible second-order correlator diagram is equivalent to 
inserting the series for the correlator, giving rise to the anticipated diagram as well as others. 



b. Construction with a Propagator Correction 



■■0(--> - a--) 



o- 



-C-\- 



(62) 



Inserting the exact propagator and its expansion obtains the desired term, (W9). 
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c. Construction with a Vertex Correction 

- ---O--- + --<1>(-^"^ ^ (63) 

As with the others, the inclusion of the vertex expansion in this case produces (W15) plus others at higher orders. 

In the above examples, the leftmost propagator and vertex have remained unrenormalized. This procedure can be 
applied to the Class-B and Class-A^ diagrams in Fig. [U generating the respective irreducible diagrams. Replacing 
the appropriate bare diagrams with the renormalized quantities in the irreducible diagrams and then collecting these 
together with the Class-A^^ terms expressed in Eq. ([57)) . leads to an equation for the exact correlator, 



+ 16-- 



+ 16-- 




Wc note the distinction between Wyld's original exact correlator equation (Fig. 5 in [3), 




where the first term on the RHS contains both A/, and A^j diagrams. It can easily be checked that both Eqs. ([64]) 
and (|65p lead to the same perturbation expansion. Further remarks about the exact correlator will be made in a later 
section; we now continue with the propagator. 

The approach Wyld used to determine an equation for the exact correlator could not be used for the exact propa- 
gator. Using both Classes A^ and A^j, the result is a primitive expansion for the propagator with some of the terms 
redundantly generated, specifically diagrams (Wp3), (Wp6), and (Wp7) in Eq. ((57)) . Wyld recognized this problem 
and circumvented this by using a Dyson equation for the propagator, expressing his arguments only mathematically 
and without a diagram equation for the exact propagator. The result was to introduce modified vertex functions 
and use the Ward-Takashi identities to relate these to the propagators, as in QFT [i^. Wyld's final complete set 
of equations has a diagrammatic expansion for the exact correlator, two diagram series for the exact and modified 
vertex functions, and the Ward-Takahashi identities. 

However, it was argued by Lee that the method of using the Dyson equation and Ward-Takahashi identities cannot 
be applied to the full 3-dimensional NSE in a manner similar to Wyld's scalar model 4J]. Lee had adapted Wyld's 
method to magnetohydrodynamic turbulence and found the same problem but introduced the following equation for 
the exact propagator: 



= = = = + 0{\^) (66) 

The left-most propagator and vertex have both been left bare. The former clears the redundant generation of 
propagator diagram (Wp3), while the latter does the same for diagrams (Wp6) and (Wp7). This equation for the 
propagator does correct the redundancies, however, the asymmetry introduced by Lee at second-order also does not 
generate Wyld's diagram (WIO) at fourth-order, hence his inclusion of a fourth-order irreducible propagator term not 
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found in Wyld. Lcc included this term in the exact propagator, 

16 :i\? = = = = + 0(A^ 



(67) 



Lee introduced this as an ad hoc fix, however we observe here that the above equation can be derived by adopting 
a scheme where the leftmost propagator and vertex of the exact propagator remain unrcnormalized, similar to what 
we did for the correlator expansion. We have now a single procedure for deriving both the correlator and propagator 
expansions. 

This was the motivation for using Eq. (|5ip at the start of this Subsection as our formulation of Wyld's approach. 
Finding the irreducible propagator diagrams from the set of Class-A^ diagrams while maintaining that the left- 
most propagator and vertex remain unrcnormalized results in Eq. (j67p . Our justification for doing so follows one 
of Kraichnan's early derivations of the propagator [l7| . His approach was to include an additional, infinitesimal 
perturbation, 7yQ.(k, i), to the equation for the velocity field (c.f. Eq. (|10p ). 

t t 

S„(k,i) = Jdt'RZ'i^-: t,t')U,{k,t') + Jdt'RZik; i,<')^a'(k,i') 

— oo — oo 



+ ^ dt'R'°l,{k; t, t')M^^Vp^{k) d'i u^j, t')u^{k - j, t'), (68) 



and then to derive the propagator as a response function to the perturbation, 
RU^;t,t')S{k-k') EE 



RZik; t,mk-k') + XRZ{k; t,t')ACZ{k) jd^\ ^-1-^ (^i,^(j, t')u^(k - j, t') 



(69) 



This response function would then be averaged in the same manner as the correlation function to give the renormalized 
response function or propagator, 



R^^{k;t,t') = (^R^^{k;t,t')y (70) 
Note that the averaging would leave the bare propagator and the vertex unchanged, 

i?„.(k;t,i')^(k-k') = i?2(k; t,t')5{k-\^) + \RZ{\^: t,s)KC'^^{k) Jd^i /^ ^^J^^ (^^^(j, g)£t^(k - j, g)^ ^ ; (71) 

this analytic form of the propagator is consistent with our current version. In particular observe that Eq. ([7T|) is 
consistent with our procedure of keeping the leftmost propagator and vertex functions as bare. 

Lastly, we consider the equation for the renormalized vertex as a function of renormalized quantities. Since there is 
no leftmost vertex or propagator, the exact vertex expansion from Wyld is still valid to 0{X^) and can be used here, 



+ 4 



+ 4 



(72) 



In summary, the exact correlator, given by Eq. (j64p along with the renormalized propagator and vertex functions, 
Eqs. (|57)) and ([7^ . respectively, constitute a closed set of integral equations that describe the time evolution of 
the two-point correlation function. We will refer to this set as the "Improved Wyld-Lee Renormalized Perturbation 
Theory". It contains the same information about the correlator expansion up to fourth-order. Fig. [1] and describes 
turbulence insofar as the NSE can be treated perturbatively, with the Gaussian statistics assumed in the external 
forcing introduced to facilitate the closure. The above equations arc still, in principle, infinite series, however the 
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rcsummation contains the detail of a greater number of terms and allows a truncation that retains more of this 
information. Truncating after second-order and noting that the bare propagator is equivalent to {dt + z/k^)~^, the 
equation above and the equation for the propagator Eq. (|67l) give the result of the DIA, Eqs. ([7|) and ([8|). 

Thus we have developed a single unified procedure for deriving both the correlator and propagator expansions in 
the Wyld formalism. The correlator equation from Wyld's original approach, Eq. (|65|). is correct, but his approach 
to calculating the propagator was flawed, as it led to the redundant generation of diagrams as observed by Lee [31 . 
Our procedure of keeping the leftmost propagator and vertex as bare recovers Wyld's original and correct correlator 
expansion but at the same time leads to a correct equation for the propagator expansion. This is our simple single 
procedure for obtaining the Wyld theory. 

As such, we have established here a well-defined procedure for renormalization. Summarizing this procedure in 
terms of the equation for the correlation, we symbolically write it as, 

C„{uu) = M'-°\uuu). (73) 

Note that the linear operator, whose inverse is the bare propagator, and the vertex function are outside of the corre- 
lations. In perturbation theory the 3-velocity correlation function above would be a function of the bare correlation, 
vertex function and response function, {uuu) = /[C'°', M'°\ and our renormalization procedure is to replace 
these by the exact respective functions /[C'"', M'"', i?'"'] f[C,M,R]. In other words our procedure keeps the 
left-most propagator and vertex out of the resummation. 



III. THE MARTIN-SIGGIA-ROSE FORMALISM 



The MSR Formalism Q is by now a well-known analytic procedure that can be used to calculate the "statistical 
dynamics of classical systems". The formalism establishes an operator theory where the observables are defined 
as Heisenberg operators. This permits a non-perturbativc treatment akin to the Schwinger formalism (46l - l48j in 
quantum field theory, which formally closes the statistical moment hierarchy. The operator formalism introduces an 
adjoint operator, which can be used in the construction of nontrivial commutation relations leading to correlation and 
response functions [l6j . Employing the Schwinger formalism for statistical closure involves the use of a generating or 
characteristic functional. An alternative to the construction of such operators comes from path integrals [iollsoj. 



Notable sources providing detailed information on MSR are works by Rose |51|, Phythian j52| - l54l |. Andersen [55[ 
and Krommes [l^. As demonstrated in their original paper, the formalism is applicable to the turbulence problem 
which has inspired further work in the analysis of least-action principles [sgI [57| and gauge symmetries [58l - l6G | in the 
study of turbulence. 

In the previous part of the paper we reformulated Wyld-Lee, in order to eliminate the double-counting problem in 
a more natural way. Then we presented a rather pedagogic clarification of their theory, particularly of the way in 
which the use of diagrams led to renormalization. Now, in this section, we aim to make a comparable elucidation of 
MSR's theory. Of course in their case there is no procedural error to rectify. Yet, a pedagogic exposition of their use 
of diagrams in obtaining renormalization is of value and, when taken in conjunction with our exposition of Improved 
Wyld-Lee theory, allows the relationship between the two formalisms to be clearly seen. In this way we are able to 
show that the two formalisms are fully equivalent. 



A. Setting up the Formalism 



It has been seen already that the velocity field Wq(x, t) is a fundamental observable in fiuid dynamics. The MSR 
formalism extends the common notion of it to that of a classical statistical operator [s^, ■ In the language of 
quantum field theory (QFT), it is similar to a Heisenberg operator, in that it is time dependent (45j . This is an 
essential first step in establishing the formalism. 



1. Dynamical Equations 



A generalized equation of motion for a generalized field variable, = UQ,j(xi,ti), is introduced, 

at,V(l) = U{1) + ;7(12)V(2) + t/(123)V(2)V(3). 

It is argued that this equation can accommodate many dynamical systems, and in principle, can be generalized to 
higher orders of interaction. We are free with the choice of variables and so to maintain consistency with the previous 
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part of this paper, we continue to work in spectral space, 

dtUaiKt) = ?7a(k;t) + [/„^(k,j;t,iO"/3(j,<0 + '^a;37(k,j,l;i,i',i'')i~i/j(j,<')^h(l>^'')- (74) 
The quantities defined within, 

Ua{k;t) o 0-point potential/external force (75) 

Uap(k,y,t,t') o 1-point potential (76) 

Uap-y(k.,j,l;t,t',t") o 2-point potential, (77) 

are the generalized interaction potential functions, which may depend on time and be random. Integration of re- 
peated arguments and summation of indices is implied. Using (5-functions for wave-vector- and time-arguments, and 
Kronecker-^s where needed, we can specify these potentials to give the spectral NSE: 

C/„(k;i) = .fa(k,0 

=^F^{k,t), (78) 
Uap{Ky,t,t')u^{lt) - -u Jd^i Jdt'S^pS{k - i)S{t ~ t')fu^{yt) 

^Dc.f!{k,y,t,t') = -u[d^} [dt'SapS{k-i)6it~t')f , (79) 



U^p^{k,yl;t,t',t")uf,{yt')u^{l,t") = AAf^^^(k) Jd^^ J d^J dt'J dt" x 

<5(k - j - l)S{t - t')S{t - t")up{y t')u^{\, t") . 
^ NCp^{k,]M ,t") ^ AAC^(k) jd']jdH jdt'jdt"5{k-i- mt - t')5{t - t") . (80) 

We see that the 0-point potential takes the role of the external force, which is uncoupled to the velocity field; the 
1-point potential comprises the linear terms of the equation, namely the dissipative term; and the 2-point potential 
now becomes the convection term. We have re-labelled the potentials to help in the comparison. 

In keeping the goal of an analogous formalism to QFT, an adjoint operator is introduced by way of a canonical 
commutation relation, 

[u„(k, t), ?Lp{k', t)] = Sa.p S{k - k'), (81) 

thus defining the operator as a functional derivative, 

Sa(k,t) = (82) 

OUa[k, t) 

In the path-integral formalism, the Fourier conjugate of the adjoint field occurs naturally in the treatment of the 
delta- functional (see for example Jensen [i^ and Krommes [31). An equation of motion for the adjoint field, 

-atS„(k,t) = 7^^„(k,j;^,^0"Mj,^0 + ^4°-;a(jJ,k;^^^^^)S^(j,^')%(l,^''), (83) 

may be constructed using Eq. ([7^ . Eq. ([ST|) and the canonical relations, 

dtua{k,t)^ — =— — - , -dtua k, t) = ; 84) 

Sua{k,t) Sua{k,t) 

the generalized Hamiltonian functional, 'H[u, u, t], is postulated and then determined from the first of these relations, 
see for example [5l|, [H, [13, [Ml ■ 

We can construct a single dynamical equation that includes Eq. ([74| and Eq. ([83|) . Both fields are collected 
together in what Eyink calls a 'doublet field' [57| and Krommes calls an 'extended field vector' [l6j. 



U{a,k;t) 



Wa(k,t) 
Ua{k,t) 



(85) 



For these 'doublet' quantities, a curly script is used and their doublet-vector indices are given in SMALLCAPS font 
where the indices are either + or — ; for example U^^+{a, k; t) = Uaik, t). We can also establish a commutator for the 
doublet field vector, 

U,{a,k-t)M.W,^-A =i[a,]J^p5{k-\^), (86) 
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where the (Pauh) matrix is 



-i 

1 



(87) 



As argument labels will soon increase, we continue with our reduced notation where all arguments are combined into 
subscript wave-vectors, WA(a,k;t) — > Wu; this is effectively the same as the notation used in MSR. 

An equation of motion for the doublet field vector is then simply constructed from the dynamical equations of u 
and u. 



Once again, the curly-script used for the so-called 'doublet-field' potentials distinguishes them from their earlier 
counterparts. The above equation combines the system of equations, Eq. (|74p and Eq. (j83p . into a single equation 
for a single field variable; it is to this equation that we apply our statistics. 



2. Statistics 



A generating functional is introduced and used to create all statistical quantities. 



cxp 



U,>,{a, k; t)rj,^{a, k; t)dt\ 



(89) 



where {• • -It denotes time-ordering. Here r}f^{a,'k,t) plays the role of the source term that is standard to these 
techniques. It is, in effect, a perturbation to the 0-point potential (c.f. Eq. ([55)) ). 

Using the generating functional, one can find the statistical moments or cumulants as needed through functional 
differentiation of the generating functional with respect to the source term. In practice, the cumulants are obtained 
by functionally differentiating the logarithm of the generating functional, returning what are called the connected 
Green's functions [i^. As an example, the first- and second-order cumulants are produced respectively via 



G^{a\ k; t) 



or]A(ct, k, t) 
{2} 



(90) 



(5r/A(a,k,i)(5r/B(/3,k',t') 
^ G3(/3;k';0- 



5ri^(a, k, t) 



(91) 



In the reduced notation, G^ia; k; t) — > Gk and Gab{<^, P', k, k'; t, t') G^ j^.'. Note that the correlator and propagator 
functions are contained within the second-order cumulant of the extended field vector, 



G(a,/3,k,k';t,t') 



?)=0 



(u„(k, t)u^j(k', t')) (ft„(k, 0"/j(k', t')) 
(S,(k,t)u^(k',t')) 



Cc,^(k,k';t,i') i?,0(k,k';t,t') 
i?/3a(k',k;t',t) 



(92) 



where we recall that Mc(k, <) is defined in 

The interest is in obtaining a dynamical equation for a particular statistical quantity, which for example can be the 
second-order correlation function of two velocity fields of a turbulent fluid. Using Eq. (|88| , Eq. ([90]) , and Eq. (|9ip , 
one can construct an equation of motion for the mean field. 



iCTj^tGk = Jic + ?yk + ^kjGj 



1 



-M 



k,j,i( '^ja + '^jd 



(93) 
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A second-order cumulant is present on the RHS in this equation for the first-order cumulant on the LHS, thus 
generating a statistical hierarchy and so the closure problem. Differentiating Eq. (j93|) by rji^r gives 

- *a,atGk,k' = 4,k' + ^^kjGj^k' + X^j,! (CjGi.k' + . (94) 

In this case, the problem of closure occurs with the last term where a third-order cumulant is introduced, 

^ = Gk,ja. (95) 

A method is needed to proceed further without the introduction of ad hoc hypotheses to link various moments or 
cumulants. 



3. Closure 



The problem of closure can now be addressed. The method employed by MSR is the Schwinger-Dyson formalism, 
which has been used to deal with the closure problem in QFT 45| . The authors sought a way to transfer this quantum 
statistical formalism to classical physics. By way of a functional Legendre transform, 

C[G,^{a;k;t)] = In Z[T],^{a;k;t)] - G,^>{a';k']t')Tj,^,{a';k';t') (96) 

a closure can be found through the introduction of vertex functions which can be related to cumulants of order-7i 
through n-order functional derivatives of the above equation with respect to Gj [l6j . As we wish to find a closure for 
the third-order cumulant, it is straightforward to calculate the three-point vertex function, 

A^ABc(a, P, 7; k, J, 1; < , t = j-^, — . ,nc^ .0 • .^xn 1 1 ^//^ 

tiGp^{a\ k; t)6G^{p\\, t')5Gc{i; 1; i ) 

^ ^G^c{P,TM',t")Y\ (98) 



5Gj,{a,k-t) 

which is symmetric with respect to the indices A, B and C. We also note this is different from usual QFT notation, in 
that we use M to represent this function instead of the more common F, which is used in both Wyld and MSR; the 
reason for doing this is for its connection to the renormalized vertex function introduced earlier. Similarly, the last 
term in (j94p may be rewritten to contain the three-point vertex 

= ^-^kliCjj'Gu'A^mjM'G^.k'- (99) 

Another function that is crucial in QFT is the self-energy function • In particle physics self-energy is responsible 
for attributing a particle with a 'dressed' or renormalized mass, which is an observable quantity. In the present 
calculation, the self-energy is defined using the three-point vertex, 

EAA'(a,a';k, k';t,t') 

= ImZ:H P, 7; k, j, 1; i, r, s)GeB' (/?, /?'; j, j'; r, r')Gec' (7, Yi 1, 1'; s')A^a'b'c' («', P' , 7'; k', j', 1'; t', /, s'), (100) 
or in reduced notation, 

Sk,k' = ^A^L^iiGjj'Gu'A^k'jM'. (101) 
The dynamical equation for the second-order cumulant, Eq. (|94p. can be rewritten with the self-energy term as, 

- *^7,5tGk,k' - 5k,k' + ^^kjGj^k' + A^k,j,iGjG,,k' + SkjGj^k'- (102) 
The inclusion of the self-energy leads to the establishment of the well-known Dyson equation [g^ 

[<k'] "Gj,k' = '^k.k' + Sk.jGj,k'. (103) 
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The Dyson equation is an equation of motion for the second-order cumulant, which is directly obtained from Eq. 
(jlOip . It also establishes the inverse bare second-order cumulant, defined by 

[G^M'] " = -^<^2dtSuM' - ^?k.k' - A^l^'ik'Cj. (104) 

Written differently, 

[Gk.k']"'= K,'k']"'-Sk,k', (105) 
this equation gives similar (Dyson) equations relating the bare and exact propagators, 

M-'=[<y-[Sk,k']±^, (106) 
which are the off-diagonal components in [G^ k'] this equation can be re-arranged to give 

fik.k' = <'k. + RZ [SjJ'] i^^j'.k'- (107) 

One can see that the renormalization is already done in the above by the self-energy term. Furthermore, this is the 
propagator (to second order) established earlier, 

i?„0(k;t,<') =<i,(k; t,t') 

+ ,(k; t, t") (^M^.s^(y)Rss- (k - j; t" , s)C^y (j; t", s) Af^,^,^, (-k)^ Rp'p{-^\ s, t'). (108) 

Returning to the vertex function, by differentiating Eq. (jlOSp and using Eq. (j98p it can be rewritten in terms of the 
self-energy and itself: 

( ^ ^kj 

'k'j' 



■^kj.i = M'Z,+ \-^^ )Gk',k"Gj'j"Mk"j'M. (109) 



This completes our brief exposition of the MSR formalism. We now have a set of equations that close the statistical 
hierarchy: 

- ia.dtG^ = J-k + Pk.yGj + \m'Z ,{G;,. + GjGz) , (110) 
-za.atGk,k' = '^k.k' + ^?k,jGj,k' + >(i:?j,iGjG,,k' + Sk,jGj,k', (HI) 
Sk,k' = ^A^LliGj,.Gu'-Mk'jM' (112) 

■^k,j,i = M^l, + (;^^)Gk',k"Gj',j"A^k"j'M. (113) 
Going one step further to incorporate the self-energy term, 

- la.dtG^ = J-k + Pk,yGj + iAlL",j,i(Gj,z + GjG,) , (114) 

-z(7.atGk,k' = ^.k' + 2?kjGj,k' + -Mkj,iGjGi,k' + (^Xk,'„,„G„,„'G„,„.A^j,„,,.n,)Gj^k', (115) 

A^kj.i = M^^, + (^-^A-(^-M'^]^^^G^^^.G^^^.^^^^^ (116) 

we find that these equations require iterations which can also extend to infinite orders. Truncations of A^kj,i are then 
necessary to deliver practical and computable results. The first-order truncation of the vertex function is 

-^k,j,i « A^l^'ii, (117) 

and while it may not be immediately evident, this truncation gives the DIA. The next order of approximation, 

A^kj.i « Mti, + M^,y,G^,^,M^P},^,G^^^,G^y,M^,.„^, (118) 

similarly generates fourth-order terms given in Wyld, as will also be seen in the following section. 
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B. The Diagrammatic Representation of MSR 

The formalism was originally demonstrated with diagrams related to the turbulence problem, and is by the author's 
account comparable to that of Wyld. Wc now work directly from their diagrammatic interpretation and make a direct 
comparison to Wyld. 

The correlator and propagators are obtained from the second-order cumulant tensor, 



G{a,l3;k,k';t,t') 



))=0 



Ca0(k,k';t,t') i?,^(k,k';i,i') 
i?^a(k',k;t',t) 



We define a diagram for this object just as with the correlator and other diagrams, 

k,t k'.t' 



B,I3 



O G,.B(a,/3;k,k';t,t') 



G 



k.k'- 



(119) 



(120) 



Using the same notation for exact correlators and propagators from the Wyld analysis, Eq. (jll9p can be transcribed 
diagrammatically as 



(121) 

The external force introduced as T would not survive in the second-order cumulant equation, Eq. ()94|) . This term 
is needed for stationary turbulence, as was seen in Wyld. MSR avoids this difficulty by introducing the forcing as 
potential. In the case of turbulence considered here, MSR considered the mean field to be zero, hence Gk = 0. Noting 
these points, Eq. (|104p is rewritten as 



[^k,k'] ' = -«cr29t(5k,k' - 2?k,k' - -^k,k'- 



(122) 



For the forcing potential 7^k,k' the only nonzero component is 
Using Eq. (jlOSp . an equation for Gk k' can be constructed. 



= (/„(k,i)/a'(k',i'))- 



GkM' — ^kj[Gjj,] Gj' k' — G'kjSjj'Gj' k' 



(123) 



From the structure of these matrices, using the Dyson equation for the propagator given in Eq. (jl07p along with 
-Rk,a[-Ra,k'] ^ = (^k.k' [iBli One Can construct a graphical interpretation for the second-order cumulant equation. 



(g, 






(124) 



Equating Eq. ()121|) and Eq. (|124|) . one can obtain a diagrammatic equation for the exact correlator 



= ===(8)=== 



==== s 




where = [^k,k']__ and 



{f{k,t)<E>f{k'X)) 



(125) 



(126) 



is the external force contribution, which in the case of NSE turbulence is the correlation of two random forces. The 
mathematical formula for this is then 



Gk,k' =i?k,j(f(j,i)®f(j',i'))%,k' +i?k,j[Sj,j.]__i?j.,k' 



(127) 
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The graphical equation for the exact propagator can hkewise be extracted from Eq. (|124|) . 



+ 




(128) 



with its analytic counterpart Eq. (|107p . Note at this point that the MSR procedure has built into it that the leftmost 
response function in Eq. (|128p is unrcnormalized and from Eq. (|10ip and Eq. (|128p that the leftmost vertex function 
is bare. Recall in the Wyld formulation Lee had introduced these corrections ad hoc, and we showed in Eq. (|5ip how 
this can be properly accounted for at the outset. With this modification to the Wyld procedure, the Wyld and MSR 
renormalized perturbation approaches arc equivalent. We will show explicitly that this is true up to fourth-order in 
the perturbation expansion. 

We now express the self-energy and vertex equations in diagrams. Using the earlier graphical notation for the bare 
and exact vertices from miBl the self-energy, Eq. (|10ip . can be presented as 



^k,k' 



k,j,k-j 




k',j',k'-j' 



k-j,k'-j' 



(129) 



Similarly, the vertex diagram, Eq. (|109p . is 



(130) 



the quantity is used here in place of (SEkj/(5Gij' j'. This equation can be substituted into the vertex term of Eq. 
(|T29| . giving 




1 
2 




G 



G 




M 



G 



where wave-vector labels have been suppressed; the corresponding equation to this diagram is 



SGm.n 

We introduce a diagram for the four-point term I4 for convenience; it can be written as 

G 



(131) 



(132) 







G _ 

7W<°' M 



G _ 

7W<°' M 




(133) 



with the corresponding equation 

h 



SGm.n 

5 (\ 



^-^j',g,h<^g,g'^h,h'-^i',g',h' 



(^Gm.n \2 
1 
2 

+ iA1<°;^_^G,,g,Gh,h'Tp^(A^i',g'.h 



-^j':m,hGh,h'A^i',n,h' + |-^^":g,™G,,g.x,^g,„ 



(134) 



24 



The term = ^A1i'.g',h'/^^m.n has also been temporarily introduced. As the next steps are intermediate, we will 
briefly suppress labels. Inserting Eq. (jl33p into Eq. (|13ip results in an expansion, 





(135) 



This can be further written with the bare vertices as 

1 

s = 





O(A^), 



(136) 



where we make the association of A with the bare vertex term as was done in the previous section. 

Note that Al'"' is a 2 x 2 x 2 tensor, and it can be shown [igI. [Slf to be symmetric with three nonzero entries. 



(137) 



Using this and noting that only [Ga.b]__ vanishes, the self-energy is constrained to have only one zero component, 
Pa,b]^_,.- The remaining terms on the RHS of the self-energy equation Eq. (|136p are found to be 




+ 



+ 



(b) 



' if) 



(d) 



(<?) 



0(A5), 



'C (h) 



(138) 



for the off-diagonal terms, and 





+ 






O(A^), 



(139) 



are the non- vanishing diagonal components. Both of these equations have been shown only to fourth-order in 
as this is the extent of this study. 

We have obtained the above results, Eqs. (|138p and (|139p from the basic equations in MSR, and as we will see 
later in this Section, they give a perturbation expansion that agrees with Wyld up to fourth order. However they 
never actually calculated in their paper the fourth order self-energy contributions as we have in Eqs. (|138p and (|139p . 
Instead MSR recast their formalism into a simplified diagramatic approach and from this their perturbation expansion 
was developed. 

The MSR diagramatic approach began by writing the fourth order contributions to the self-energy in terms of two 
Green's function lines convoluted with vertex functions. They wrote the non-vanishing elements of the self-energy 
tensor through the diagrammatic equations 









(140) 




(141) 
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where a, /3, and 7 label the three vertex functions, which MSR claimed to be 




(142) 



(143) 



(144) 



These three terms correspond to the non- vanishing elements of A^, which are [A^]+^_, [A^]__+, and [A4] , and their 

permutations. We have displayed these vertex functions noting specifically that for a given combination of indices, 
the a-vertex function has four diagrams, the /3-vertex has two diagrams, and the 7-vertex only one. The expressions 
on the LHS in the above equations indicate how vertex diagrams attach to the self-energy terms. Permuting indices 
leaves the vertices unchanged, for example [A^]__+ = = [A^]+__, and diagrammatically this corresponds 

to in-plane rotations and reflections about a line connecting a vertex and the midpoint of an opposing edge. Such 
permutations are trivial for the a- and 7-vertex functions where any operation returns the same diagrams given in the 
above equations. The /3- vertices require these permutations to generate the necessary self-energy terms. For example, 
E:p± uses Eq. (|143p as it is shown; however, requires a 120° counterclockwise rotation, [A^]__^ -> [A^]_^„„, for 
diagrams Eq. (|139[) (' f) and (e), respectively, whereas a 120° clockwise rotation, [A^]__+ is needed for Eq. 

p39p fo') and (h), respectively. In this manner, one can confirm that all the diagrams in Eqs. p38p and (|139p are 
produced. We will examine now in close detail how the formalisms of MSR and Wyld lead to complete agreement 
at fourth order. Inserting the vertex corrections, Eqs. (|142p . (|143p . and (|144p into their respective positions in the 
self-energy diagrams, Eq. (|140p and Eq. ()14ip . gives the self-energy diagrams to fourth-order. Then, the expanded 
self-energies are inserted into Eq. (|128p to obtain from the MSR formalism the propagator diagrams to fourth-order. 




A factor p has been introduced for each vertex since this overall factor differs between MSR and Wyld. Using this 
propagator expansion, the terms in Eq. (|125p for the exact correlator may now be determined. Those terms obtained 



from the forcing function are 
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(146) 



The labels (Wn) correspond to Wyld's diagrams as seen in Fig. [TJ Labels containing an asterisk denote diagrams 
with half the weighting of their Wyld counterparts; however there are always two such diagrams and their sum gives 
the correct weighting. These diagrams are symmetric when reflected about a horizontal line and in Wyld's formalism 
such diagrams are equivalent. 
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Those diagrams representing the self-energy interaction are expressed by 




(147) 



Combining Eq. (|146|) and Eq. (|147p together gives the diagrammatic equation for the exact correlator function 
expanded to fourth-order in the bare vertex as obtained by the MSR formalism. There are 44 diagrams counted in 
these two equations, confirming a statement made in the MSR paper. Some of these 44 diagrams are repetitions of 
the same diagram. If these diagrams are then combined taking appropriate account for weight factors, they yield the 
33 diagrams in the Wyld perturbation expansion in Fig. [TJ Comparing these two equations to the diagrams given in 
Fig. [U it can be concluded that the MSR formalism gives the primitive correlator expansion of Wyld provided the 
factor p — 2. 

IV. COMPARISON 

The renormalized diagrammatic expressions for homogeneous isotropic turbulence have been presented here accord- 
ing to the formalisms of Wyld (sections '^ ^W BfHTCj) and Martin, Siggia, and Rose ( ijlll Bp . Comparisons arc made 
between the Improved Wyld-Lee formalism presented in section §11 CI to expressions for the self-energy and vertex 
diagrams derived by MSR in the previous section. Ultimately we show that both formalisms are equivalent to fourth 
order. 

To compare the diagram equations as obtained in their respective formalisms, some rearrangement of terms is 
necessary. Looking again at the propagator expression of Wyld in Eq. (j67p . this equation can be written in such a 
way as to anticipate the MSR form of the propagator with the self-energy term 




(148) 



Noting the reflection-symmetry about the horizontal line of Wyld diagrams, the fourth-order term in Eq. (|148|) can 
be written as a sum of two terms. 




(149) 
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This allows the above equation for the propagator, Eq. (I148|) . to be written using two of the MSR- vertex terms, ^A^°'^ 
and M'^^^ (see Eq. ([Ti^ and Eq. p^ ): 



© + 0(A6) 



(150) 



In all further diagrams of both formalisms, diagrams that are equivalent through this symmetry will be combined. 
The "©(A^)" has been kept explicit as the above only applies to fourth-order; the next order has not been analysed. 
Substituting the off-diagonal self-energy term, Eq. (|14ip . in place of the terms in the brackets of Eq. (|150p . we find 
that the MSR propagator given in Eq. (|128p is recovered, 




Before continuing with Wyld, it will be helpful to rewrite the diagram equation for the exact MSR-corrclator, seen in 
Eq. (|125p and reproduced here, 



= = = = S 




Substituting the propagator equation given above (from Eq. (|128p ) for the exact propagator on the left of the 
force-force correlation, (f (55 f) and the self-energy term, S_„, the correlator equation becomes 







(151) 



However, in making the comparison with Wyld, it will be better to use Eq. (|15ip . To fourth-order, the self-energy 
terms can be written explicitly as 





.4 




(152) 




= p- 



A o 



©(A'^) 



(153) 



The factor p, associated with each vertex, has been included for convenience. Note that independent of MSR, these 
terms had been derived at about the same time as them for the case of wave turbulence in [64 1. The terms with only 
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one self-energy component are given as 



and 




(154) 



(155) 



The term with both !]__ and will not be expanded. This term can be combined with Eq. (|155p to give 




(156) 



Using the expressions Eq. (|154p - ()156p . these equations can be substituted back into Eq. (|151|) . 




(157) 



The term p — 2 has been used, which will allow the comparison with Wyld. To make a comparison with this correlator, 



30 

some adjustment is required for the Wyld correlator. The equation for the correlator given in Eq. (|64p . 




(158) 



can be expanded using the renormalized vertex expansion of Eq. (|58[) and truncated to fourth-order, leaving 




(159) 



This equation can be readily compared to Ec^. (jl57p . The double- force correlation term in Wyld can be connected 
with that given in MSR (see Eq. (IT26)) ). 

iff) = (g) (160) 



The equations Eq. ()157|) and Eq. (|159|) are equivalent. It can be said that the results produced here by both 
formalisms are equivalent to fourth-order. This raises questions on the claims by MSR and Kraichnan over the incorrect 
treatment of the vertex renormalization by Wyld. Their claim was that Wyld was missing certain classes of vertex 
corrections, namely Eq. (|144p and Eq. (|143p . This is true, but Wyld's resummation procedure is also different from 
MSR. In the Wyld approach, at every order in the perturbation expansion, he adds what he calls irreducible diagrams. 
Notice in the original Wyld formulation his correlation equation (|65|) had two types of irreducible diagrams; one can 
compare to the corresponding MSR equation (|140p and see the difference. In making this comparison, importantly 
note that the leftmost vertex in the Wyld case is also renormalized, which again differs from MSR. Nevertheless, 
as we saw on comparing Wyld and MSR, that up to fourth order, they lead to equivalent perturbation expansions. 
We made a modification to the original Wyld procedure, in order to have a single procedure for treating both the 
propagator and correlator. For this, we modified the correlator equation of Wyld to Eq. in which, just like for 

the propagator, the leftmost vertex and propagator remain unrenormalized. However the basic principle adopted by 
Wyld to include certain irreducible diagrams at each order in the perturbation expansion is still adhered to in our 
Improved Wyld-Lee perturbation theory. Within this approach, again we explicitly confirmed that Eq. (j64[) agreed 
with MSR to fourth order. 
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V. CONCLUSION 



The formalisms of Wyld and Martin, Siggia, and Rose have been presented in detail is this paper. Specific attention 
has been given to the diagrammatic interpretations of both formalisms, as this is where confusion can arise when 
making comparison. 

A new derivation of the propagator and correlator for the Wyld formalism has been given, that obtains the corre- 
sponding expressions proposed by Lee. The main feature of this resummation procedure is that the leftmost propagator 
and leftmost vertex remain unrenormalized. This procedure, when added to the original Wyld formalism, we named 
here as the Improved Wyld-Lee Renormalized Perturbation Theory, which in particular comprises Eqs. ([M)) . (|67p . 
and ([72]) . Applying this method for the propagator, one directly obtains the DIA results when truncated to lowest 
nontrivial order, just as one does with the MSR formalism. 

Furthermore, equations of the Improved Wyld-Lee Formalism for the exact correlator ([M)) and propagator (p7)) 
have been shown to agree with the corresponding MSR equations (|140[) and (|141|) respectively. What we have pointed 
out here is that Wyld had a different procedure for expressing his basic equations, involving what he called irreducible 
terms. Accounting for these, we illustrated how the Wyld procedure agrees with MSR up to fourth order. We had 
included the factor p in the MSR diagrams, (jl54l) and Eq. (|155p . to account for a weighting factor. We found that 
setting p = 2 gives the proper Wyld weightings but this cannot be known a priori by the MSR procedure. However our 
corrections do not fundamentally change Wyld's approach of classifying certain irreducible diagrams and accounting 
for them in his perturbation expansion. It is this procedure that makes his approach different from MSR, and due 
to that, just a simple comparison of the vertex functions in both formalisms is insufficient in making a comparison of 
the two formalisms. It is therefore the conclusion of this work that both formalisms are equivalent. This means that 
the Wyld formalism, with the diagrammatic resummation used here, produces equations for the exact correlator and 
propagator that are the same as those obtained using the formalism of MSR. 
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